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Abstract 

We present a theoretical study of near-field heat transfer between two uniaxial anisotropic planar 
structures. We investigate how the distance and relative orientation (with respect to their optical 
axes) between the objects affect the heat flux. In particular, we show that by changing the angle 
between the optical axes it is possible in certain cases to modulate the net heat flux up to 90% at 
room temperature, and discuss possible applications of such a strong effect. 
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Since the prediction by Polder and van Hove [ij that the heat exchange between two media 
at short separations can be much higher than the blackbody hmit, numerous works have 
been carried out to investigate both theoretically and experimentally the physics involved 



Moreover, very recent experiments [6|, 3] were in good quantitative agreement with 



in this transfer. Experimentally, it was shown |2|, |3[ that the radiative heat flux increases for 
distances shorter than the thermal wavelength and can vastly exceed the black body limit 

a. 

theoretical predictions. On the theoretical side, we can highlight the studies of the heat flux 
for layered media 3; ^or photonic crystals lo|, metamaterials 11|, and porous media 12|. 
In addition, the dependence of the heat transfer on the geometry has attracted much interest 



and has been investigated in a sphere-plane geometry 
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plane surface 



15| and between two spheres or nanoparticles 



, for s 
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^heroidal particles above a 



42 of . Somewhat more applied 



studies have attempted to take advantage of the potential of the tremendous increase of the 
radiative heat flux on the nanoscale for thermal imaging of nanostructured surfaces 



Finally, the formulation of the heat flux in terms of the scattering matrix 
the way for the study of further geometries and the Landauer concept 27, 
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24|. 



26[ paves 



opens up a 



deeper understanding of the trade-off between heat transmission and the number of modes 
contributing to the heat flux. 

While considerable progress has been made over the last decades to actively manage heat 



flow carried by phonons in nanostructures 



29 



-|31|. very few attention has been paid so far 



on the control of non-contact heat exchanges at nanoscale. In 2010, a pioneer work carried 
out in this way by Otey et al. proposed a thermal rectifier based on photon tunneling 
between two thermally dependent polar materials separated by a vacuum gap. The efficiency 
of the thermal rectification was found to be about 40%. More recently, a device made with 



phase change materials has been introduced by van Zwol et al. [33| to modulate heat fiux 
between two materials using an electric ac current as external power source. However, due 
to the properties of these materials, such modulator works reversibly only during a limited 
number of cycle which typically oscillate between 10^ and 10^^. Moreover, such devices work 
at two discrete levels of fiux, one for each state of the phase changing material. 

In this Letter, we investigate the near-field heat transfer between two polar/metallic 
misaligned gratings in the long wavelength limit, where they may be described by effective 
homogeneous anisotropic permittivities. We show that it is possible to get a strong heat fiux 
modulation without cycle limitation just by rotating the relative position of the grating's 
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optical axes 3J]. Our approach combines the standard stochastic electrodynamics 35|] and 



the effective medium theory 



36- 



for the gratings. 

A sketch of the geometry considered is depicted in Fig. [TJ It shows two semi-infinite host 
materials of complex permittivity €i{uj) {i = 1,2) with a one dimensional grating engraved 
on each. The relative orientation of the two gratings is arbitrary in the (x,y) plane, and 
we assume that their trenches are sufficiently deep so as to (i) render the substrate below 
those gratings irrelevant and (ii) allow us to consider the x- and z-directions as equivalent. 
Moreover, these structures are separated by a vacuum gap of thickness d and kept at two 
different temperatures Ti, T2 in local thermal equilibrium. By choosing a coordinate system 
with the z-axis perpendicular to the interfaces, we can write the permittivity in the form 
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€±. sin (pi + ey^ cos < 
(e±, - e||J sin 0i cos I 



le±, - eilJsm^i cos> 
cos^ (-^ ' ' "'"^ 



(1) 





+ e|j, sin^ (pi 
ey 

is the angle between the y- and the i-th structure's optical axes, and en, e± are 



where 

given by the effective medium expressions 



36 



eh,: 



(2) 



:i - f^) + /.eh, 

where /j and e/j. are respectively the vacuum filling factor and permittivity of the host 
material in the i-th grating. These expressions are valid for arbitrary filling factors as long 
as the grating periods Ai are much smaller than the thermal wavelength Ath = ch/ksT, 
that in our case (300 K) is about 7.68 fim. Within the near-field regime this condition for 
the validity of the homogenization is different. As discussed in more detail in Ref. 12| the 
expressions for the effective permittivity in the near-field regime are valid if Aj is smaller 
than the spacing d between the gratings. 

According to standard stochastic electrodynamics jl|,[35|, the heat flux exchanged between 
the two bodies per unit surface is given by the statistical average of the Poynting vector 
normal component 

dk 



where Q{u,T) 



r^mu,T,)-Q{u,T,)] 

Jo 



(27r) 



(3) 



1] is the mean energy of a harmonic oscillator and the 
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transmission factor T(a;, fcy) [27| can be written as [12 1 



Tr \(1 - R^R2)D^^(1 - RiRI)]D^2^ ^ n<u/c 
T{u,H.;d) = { [ J (4) 

(RJ-R2)D^2(Ri-R1)D^^^ e-^l^l'^ K>uj/c, 



Tr 



for propagating (k < u/c) and evanescent (k > w/c) modes where 7i. = ^yu'^/c^ — k^, and 
= (l-RiR2e2i>'^)-i. The reflection matrix R^ of the i-th structure is a 2x2 matrix in the 
polarization representation. Its four elements Rm with k, I G {s, p} for the scattering of s- or 
p-polarized plane waves into s- or p-polarized for the considered structures are determined 



with the method presented in Ref. 40 1. 



Now, we discuss the numerical results obtained with the above expressions for two Au 
gratings. In Fig. |2] (a) we show A (5,) = [(5,(0 = 0°)) - (5,(0 = 9O°))]/(5,(0 = 0°)) over the 
filling factor for three different distances d = 100 nm, 500 nm and 1000 nm. Here, </> = 0i — 02 
is the twisting angle between the optical axes of the gratings. One can observe that the 
difference in the heat flux between both configurations, i.e., two parallel or perpendicular 
gratings, is larger than 70% for all three distances and for a wide range of filling factors 
ranging from / = 0.1 to 0.9 and can even reach about 90% for d = 100 nm. In Fig. |2] (b) 
the dependence of (5,(0)) on the twisting angle for a fixed filling factor of / = 0.3 shows 
that for all considered distances the heat flux is very sensitive to a twisting of the gratings. 
The flux drops at least by 50% for relatively small twisting angles of = 30°. We note that 
(5,(0 = 0)) is about 12, 3.8, and 2.5 times the black body value ((5bb) ~ 459.3W/m2) for 
d = 100 nm, 500 nm, and 1000 nm and / = 0.3. 

The observed sensitivity of the heat flux with respect to the twisting angle can be 
understood by the following consideration. When we have two parallel Au gratings, both 
still have a metallic response for a plane wave with polarization parallel to the grating 
structure. On the other hand, for a plane wave perpendicularly polarized to the gratings we 



have, according to Eq. ([2]), a dielectric-like response dictated by e± [4l|. For plane waves with 
polarizations between these two cases we have a mixture of metallic and dielectric response, 
but the crucial point is that we have always a symmetric situation which favors a larger 



heat flux |42|. By twisting the two Au gratings we break this symmetry. In particular, 
it is impossible to have a simultaneous metallic response of both gratings for any fixed 
polarization. 

Finally, we will focus on the numerical results for two dielectric gratings made of SiC. 



In Fig. |3] (a) we show the same plots of A{Sz) as in Fig. |2] (a) but for two SiC gratings. 
Obviously, in this case the difference in the heat flux in the parallel and perpendicular 
configuration is smaller than for the two gold gratings and it varies much more with respect 
to the filling factor. The heat flux is in this case also less sensitive with respect to the 
twisting angle as shown in Fig. [3] (b). We note that {Sz{4' = 0°)) is about 22.9, 4.2, and 
2.8 times the black body value for respectively d = 100 nm, 500 nm, and 1000 nm, with 
/ = 0.3. While the drop in the heat flux when twisting the gratings can still be ascribed to 
the breaking of the symmetry, the underlying physical mechanisms are more involved, since 
in contrast to Au gratings for SiC gratings coupled surface modes and frustrated modes 
determine the heat flux at T = 300 K. By changing the filling factor the mode structure 
of these surface and frustrated modes is changed and we find that new surface modes and 
bands of frustrated modes appear as was also found for porous media |l^. By twisting the 
gratings the coupling between these modes becomes less efficient resulting in a smaller heat 
flux. 

In conclusion, we have theoretically shown that the near-field heat flux exchanged between 
two parallel grating structures can be modulated up to 90% by acting on the relative position 
of the optical axes of the gratings. We have also demonstrated that the flux magnitude is very 
sensitive to the twisting angle, even for a low filling factors. This allows for manipulating 
the heat flux at nanoscale which can be usefuU for thermal management in micro/nano 
electromechanical devices. On the other hand, an efficient active control of transmission 
properties by mechanically driven grating structures might be impractical. However, one can 
expect a similar effect for materials such as liquid crystals or metal ferromagnetic structures 
for which the optical axis can be easily controled by applying external fields. 
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Figure 1: Two gratings separated by a distance d, and relatively twisted by an angle (p. 
(a) ■ 
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Figure 2: The heat flux between two Au gratings, (a) A{Sz){f) is defined as the difference in the 
heat flux between two parallel and two perpendicular gratings, normalized by the former (see also 
text), (b) normalized by the flux {Sz){0°) when the gratings are aligned. The angle (p measures 
the relative twisting between the gratings, and the filling factor is fixed at / = 0.3 





Figure 3: The same plots as in Fig. [2l but for SiC instead of Au. 
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